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APPLICATIONS  OF  ALTERNATIVE  PROBLEMS 


Jack  K.  Hale 

1.  Introduction.  Many  problems  in  analysis  and  applied  mathematics 
can  be  reduced  to  the  solution  of  equations  in  a  function  space  or 
functional  equations.  The  equations  often  arise  from  the  desire  to 
obtain  solutions  of  ordinary  or  partial  differential  equations  with 
subsidiary  conditions  -  periodicity  or  more  general  boundary  con¬ 
ditions,  specified  asymptotic  behavior,  analyticity  conditions,  etc. 
Many  of  the  problems  involve  a  linear  operator  and  a  nonlinear  opera¬ 
tor  which  is  small  when  some  parameter  is  small.  If  the  linear 
operator  has  an  inverse,  conceptually  there  are  no  difficulties  in 
obtaining  approximate  solutions  although  practically  there  may  be 
many  difficulties.  At  least  one  can  draw  on  all  existing  fixed 
point  techniques .  When  the  linear  operator  has  elements  in  its  null 
space,  then  new  concepts  must  be  introduced  in  order  to  proceed. 
Generally,  one  relies  on  some  type  of  "Fredholm"  alternative  to 
determine  conditions  on  the  elements  of  the  null  space  of  the  linear 
operator  in  order  to  obtain  solutions  of  the  nonlinear  problem.  For 
equations  with  a  small  parameter  in  the  nonlinearity,  one  often  ex¬ 
pands  everything  in  a  power  series  in  the  parameter  and  equates 
coefficients  in  the  spirit  of  the  method  of  Poincar/  for  obtaining 
periodic  solutions  of  ordinary  differential  equations .  On  the  other 
'  and,  experience  over  the  ldst  twenty  years  has  shown  that  it  is 
advantageous  to  look  at  these  latter  questions  in  a  more  general  way 
in  order  to  understand  the  underlying  structure.  In  fact,  proceeding 
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in  this  fashion  leads  not  only  to  new  results  hut  also  to  other 

iterative  methods  for  obtaining  approximations  to  the  solutions. 

It  is  purpose  of  these  lectures  to  present  a  general  approach 

to  the  solution  of  functional  equations  and  to  indicate  how  to  obtain 

the  appropriate  functional  equations  for  a  variety  of  applications . 

Let  X,  Z  be  Banach  spaces  and  B!  X  -*  Z  be  a  mapping  from 

same  subset  of  X  into  a  subset  of  Z.  The  range,  domain,  and  null 

space  of  B  will  be  denoted  by  j^(B),  ^(B),  J}(B),  respectively. 

A  projection  R  in  a  Banach  space  X  is  a  continuous  linear  mapping 

2 

taving  X  into  X  such  that  R  =  R.  If  R  is  a  projection  in  X, 

the  range  of  R  will  be  denoted  by  X^  and  the  symbol  X^  will 

always  denote  a  subspace  of  X  which  is  obtained  through  a  projection 

operator  R.  The  identity  operator  will  be  denoted  by  I.  If  R  is 

a  projection  in  X,  then  the  space  X  is  a  direct  sum  of  X„  and 

X  and  x  e  X  can  be  written  in  a  unique  manner  as  x  =  x +  x_ 

_L -K  K  X  ~K 

"lth  *R  6  XE’  *I-E  €  *I-R- 

Let  A.'  ^(A)  C  X  -»  Z  be  a  linear  operator  defined  in  the 
subspace  and  R!  X  -»  Z  be  a  linear  or  nonlinear  map.  The 

basic  problem  is  to  find  an  x  e  Q>(A)  such  that 


Ax  =  Nx. 


(1.1) 


In  later  sections,  we  show  how  this  type  of  equation  arises 
in  the  theory  of  bifurcation  in  integral  equations,  partial  differ¬ 
ential  equations  and  ordinary  differential  equations,  the  theory  of 
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differential  equations  with  singular  points,  the  theory  of  boundary 
value  problems  and  asymptotics.  For  the  next  few  pages  we  will  be 
concerned  primarily  with  methods  of  solution  for  (1.1).  The  spirit 
of  the  general  presentation  in  Section  2  has  its  origin  in  a  paper 
of  Cesari  [4] .  Personal  conversations  of  the  author  with 
Stephen  A.  Williams  led  to  the  particular  formulation  given  below. 
The  author  is  indebted  to  Robert  Glassey  and  Orlando  Lopez  for 
assistance  in  the  preparation  of  the  notes. 


2.  General  Theory.  As  we  have  pointed  out,  our  main  purpose  is  to 


solve  I  Ax  =  Nx,  where  A  is  a  linear  operator  defined  on  some  sub¬ 
space  3>(k)  of  a  Banach  space  X  with  values  on  Z  and  N  is  a 
linear  or  nonlinear  operator. 

Let  us  assume  the  following'. 

(Hi)  there  exist  projections  U  and  E  such  that 
9d(A)  =  Xy  and  &{k)  =  Z^ 

Clearly  this  assumption  implies  A  is  one-to-one  on  X^.  y  fl  &(k)  . 
Furthermore,  for  each  x  e  3>{k) ,  x  =  Xy  +  x^  we  have  Ax  =  Ax^-  y 
because  AU  =  0  and  so  the  image  of  x  ^  H  &(k)  under  A  is 
jS?(A).  This  means  there  exists  a  linear  operator  M'.  &(k)  -»X 
such  that  AM  =  I  on  &(k),  M  is  onto  Xj.  D  &(k)  and  MA  =  I 
on  Xj  y  fl  3f{k).  Then,  for  each  x  e  &(k) ,  we  can  write'. 

MAx  =  MAXy  +  MAXj  y  =  xj  y  =  (i-U)x  V  x  e  &>(k) 

and  so  MA  =  I  -  U  on  Q>(k) . 

Obviously,  the  equation  Ax  -  Wx  =  0  is  equivalent  to  the 
system*.  (l-E)(A-N)x  =  0  and  E(A-N)x  =  0.  Noticing  that  (i-E)A  =  0, 
EA=A,  MA  =  I  -  U  we  obtain  the  following! 

Lemma  2.1.  If  (Hi)  is  satisfied,  then  there  exists  a  linear  operator 
Ml  Zg  -»  3>(k)  such  that  AM  =  I  on  Z^,,  MA  =  I  -  U  on  &(k)  and  the 


equation 
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Ax  =  Nx 


is  equivalent  to  the  system 

a)  x  =  Ux  +  MENx 

b)  (i-E)Nx  =  0. 


(2.1) 


(2.2) 


Any  method  for  the  determination  of  a  solution  of  (2.1)  must 
take  into  account  Lemma  2.1.  In  fact,  equation  (2.2b)  says  that  Nx 
must  be  in  the  range  of  A  and,  if  Nx  is  in  j^(A),  (2.2a)  says 
that  the  solution  of  (2.1)  is  a  particular  solution  MENx  plus  an 
element  of  the  null  space  of  A. 

Lemma  2.1  is  basic  for  many  of  the  known  methods  for  solving 
problems  which  involve  a  small  nonlinearity  (say  the  nonlinearity  is 
continuous  in  a  small  parameter  and  vanishes  together  with  its  first 
derivative  for  the  parameter  equal  to  zero).  More  specifically,  one 
can  fix  an  arbitrary  element  y  e  9i(A)  and  solve  the  equation! 


x  =  y  +  MENx  (2.3) 

for  a  function  x  =  x*(y).  The  function  x*(y)  will  be  a  solution 
of  (2.1)  if  y  can  be  determined  in  such  a  way  that 


(l-E)Nx*(y)  =  0. 


(2.4) 
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These  latter  equations  are  usually  referred  to  as  the  bi¬ 
furcation  equations  or  determining  equations  for  (2.1).  After 
careful  study,  the  method  of  solution  indicated  in  (2.3),  (2.4)  will 
be  seen  to  be  the  underlying  principle  in  the  papers  of  Cesari  [1  -  3, 
11],  Hale  [1  -  4],  Perel^  [1],  Friedrichs  [1],  Cronin  [1],  Bartle  [1], 
Graves  [1],  Lewis  [1,  2],  Bass  [1,  2],  Nirenberg  [1],  Vainberg  and 
Tregonin  [1],  and  Antosiewicz  [1,  2]. 

As  we  shall  see  later,  many  different  iterative  schemes  can 
be  devised  for  the  successive  determination  of  approximations  to 
x*(y)  and  y. 

Another  method  for  attempting  to  solve  equations  (2.2)  is 
to  write  them  as 


a)  xx_v  =  MEN(xu  +  x]._u) 

b)  (i-EWXy  +  xx_v)  =  0. 


(2.5) 


For  a  fixed  y  e  x^  ^  one  can  try  to  determine  x^.  =  x^(y) 

so  that 


(i-E)N(x^(y)  +  y)  =  0  (2.6) 

and  then  determine  y  so  that 

y  =  MEN(xJ(y)  +  y).  (2.7) 
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This  latter  approach  is  very  similar  to  the  one  used  by  Rabinowitz  [1] 
and  Hall  [1]  for  a  problem  in  partial  differential  equations. 

Let  us  make  our  second  assumption! 

(H2)  there  exists  a  projection  S!  X  -> X  with 
X  C  X  0  ^(A)  and  SU  =  0. 

o  X—  U 

Of  course  the  assumption  implies  that  US  =  0  and  so  R  = 

I  -  U  -  S  is  a  projection. 


Lemma  2.2.  Suppose  (HI),  (H2)  are  satisfied,  R  =  I  -  U  -  S,  M  is 
the  right  inverse  of  A  as  in  Lemma  2.1  and  x  e  X  is  written  as 
x  =  xu+x 


V 


Then  the  equation  (2.1)  is  equivalent  to 


a)  x^  -  RMEN(xu  +  xg  +  xR)  =  0 

b)  xg  -  SMEN(xu  +  Xg  +  x^)  =  0  (2.8) 

c)  (l-E)N(xu  +  Xg  +  xR)  =  0 


Proof!  Lemma  2.1  implies  (2-1)  is  equivalent  to  (2.2)  and  SU  =  0 
implies  that  (2.2a)  is  equivalent  to  (2.3a),  (2.3b).  This  proves 
the  result. 


Lemma  2.3 .  Under  the  hypotheses  of  Lemma  2.2^  there  is  a  constant 
k  >  0  such  that  |RM|  <  k  (or  | SM]  <  k| )  if  and  only  if  |Rx|  < 
k|Ax|  (or  |Sx|  <  k|Ax|)  for  all  x  e  Q>(k) . 
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Proof!  For  any  zeZ  and  x  6  3>(A)  for  which  Ax  =  z,  we  have 

RMz  =  RMAx  =  R(l-U)x  =  Rx  because  RU  =  0  (or  SMz  =  SMAx  =  S(l-U)x 
=  Sx  because  SU  =  0).  The  result  is  now  obvious. 

Corollary  2.1.  If  (Hi)  is  satisfied,  then  A  has  a  bounded  right 
inverse  M*.  Xj  y  n  2>(A)  if  and  only  if  there  is  a  constant  k 

such  that 

|  (l-U)x|  <  k|Ax|  for  x  e  2>{A) .  (2>9) 

Furthermore  | (i-U)M]  <  k. 

Proof.  Take  S  =  0  and  apply  Lemma  2.3. 

This  corollary  shows  in  particular  that  the  computation  of 
the  right  inverse  is  unnecessary  since  a  bound  can  be  obtained  using 
only  the  known  operator  A. 

The  existence  of  a  bounded  right  inverse  for  A  can  often 
be  deduced  from  the  following  result. 

Lemma  2.4.  If  (Hi)  is  satisfied  and  A  is  a  closed  operator,  then 
A  has  a  bounded  right  inverse . 

For  a  proof,  see  Nirenberg  [1], 

Lemma  2.2  is  the  simple  observation  that  equation  (2.2a)  may 
be  written  as  two  equations  by  means  of  the  projection  operators  R,  S 
and  is  applicable  to  problems  in  which  N  may  not  necessarily  be  small. 
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By  choosing  the  operator  S  appropriately,  for  a  fixed  y  €  X^, 
z  e  one  can  attempt  to  solve  the  equation 


xR  -  RMEN  (y+z+x^)  =  0  for  an  (2.10) 

for  an  x^  =  x^(y+z)  and  then  determine  y,  z  so  that 

z  -  SMEN(y+z+x*(y+z))  =  0 

(2.11) 

(l-E)N(y+z+x*(y+z))  =  0  . 


A  special  case  of  this  procedure  can  be  traced  to  the  paper 
of  Cesari  [4].  Knowing  a  priori  that  x^(y+z)  is  small  for  y,  z 
in  some  bounded  set^  then  a  natural  first  approximation  for  the  solu¬ 
tion  of  (2.1)  is  to  let  x^  =  0  in  (2.11)  and  solve  the  resulting 
equations.  Retracing  the  steps  through  which  we  arrived  at  (2.11)^ 
one  observes  that  this  approximation  corresponds  to  a  generalized 
Galerkin  approximation.  The  general  procedure  outlined  above  can 
actually  be  used  as  a  theoretical  way  to  justify  Galerkin T s  method 
(see  Cesari  [4] ) . 

In  spite  of  the  simplicity  of  Lemma  2.2^  it  seems  to  in¬ 
clude  the  method  given  by  Bancroft ^  Hale  and  Sweet  [1]  as  far  as  the 
applications  are  concerned.  To  see  this^  we  give  a  detailed 
description  of  the  analogue  of  Lemma  2.2  in  the  cited  paper.  Suppose 
(HI)  and  (H2)  are  satisfied^  A  has  a  bounded  right  inverse  M!  Z„ 
n  0(A)  and  let  P  =  U  +  S.  Then  there  is  a  projection 
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QI  Z  -»  Z  such  that  (2.1)  is  equivalent  to  the  equations 


a)  x  -  Px  +  MQJIx 

b)  (i-Q) (A-N)x  =  0. 


(2.12) 


These  equations  are  very  interesting  because  they  have  the 
same  form  as  the  equations  (2.2)  since  (2.2b)  is  equivalent  to 
(l-E)(A-N)x  =0.  It  can  also  be  shown  that 


(I-Q)A  =  AP,  I-Q=I-E  +  Q  where 
(i-E)Q  =  Q(I-E)  =  0. 


Therefore^  if  x  =  x^  +  Xg  +  x^,  R=I-U-S=I-P, 
then  (2.12b)  is  equivalent  to 


Axg  -  QN(xu+Xg+xR)  =  0 
(l-E)W(xu+Xg+xR)  =  0. 


Since  MZ  =  X~  fl  Ql(k) }  it  follows  that  (2.12)  is  equivalent  to 

Q 

a)  xR  -  MQN(xu+xs+xR)  =  0 

b)  Xg  -  MQE(xu+Xg+xR)  =0  (2-13) 


c)  (i-.Q) (A-E) (Xg+Xg+Xg)  =  0. 
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These  equations  have  the  same  form  as  (2.8)  except  with  RM  =  MQ, 

SM  =  MQ.  On  the  other  hand,  Lemma  2.3  shows  that  the  norm  of  these 
operators  can  be  estimated  in  terms  of  the  operator  A  so  that  it 
is  unnecessary  to  discuss  the  operator  Q. 

We  now  give  some  sufficient  conditions  for  the  solvability 
of  (2.2a) ,  and  for  that  we  need  the  following  assumptions! 

(H3)  there  is  a  constant  k  >  0  such  that  |Rx|  < 
k| Ax|  for  all  x  e  3>{k) ,  where  R  =  I  -  U  -  S. 

(h4)  there  are  positive  continuous  nondecreasing 
functions  a(p),  p(p),  0  <  p  <  ~  such  that 

M  <  p(P) 

| Nx-Ny|  <a(p)|x-y|  for  ] x]  <  p,  |y|  <  p. 

For  any  positive  constants  c,  d  let 

V(c)  =  {x  e  Xj^!  |x|  <  c) 

S(d)  =  {x  e  XR!  | x |  <  d] 

Theorem  2.1.  Suppose  (HI)  -  (h4)  are  satisfied,  M  is  the  operator 
in  Lemma  2.1  and  c}  d  are  such  that 


RME|a(c+d)  <  1,  | RME| p (c+d)  <  d 


(2.14) 
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Then  for  any  x  in  V(c),  the  map 
A(-)  =  RMEN(x  +  •):  Xjj  ->X 

is  a  contraction  mapping  of  S(d)  into  itself  with  a  resulting  unique 
fixed  point  Fx.  The  function  F!  V(c)  -> is  continuous.  Further¬ 
more,  if  N  has  a  continuous  Frechet  derivative,  then  F  has  a 
continuous  Frechet  derivative. 


Proof  I  Lemma  2.3  implies  |RME|  <  lc|  E  ]  is  hounded  and  the  result 
is  now  easily  verified  using  the  contraction  principle. 


Theorem  2.2.  Suppose  the  conditions  of  Theorem  2.1  are  satisfied  and 
F  is  the  function  given  in  Theorem  2.1.  Then  the  equation  Ax  =  Nx 

•Xy)  €  V(c)  and 


has  a  solution  if  there  exist  xg,  x^  such  that  (xc 


a)  xg  -  SMENCXy+x^FCXy+Xg))  =  0 

b)  (i-E^Xy+Xg+F  (Xy+Xg))  =  0. 


(2-15) 


Conversely,  if  there  exists  an  x  such  that  Ax  =  Nx,  Jx^l  <  d, 

lxU+Xs'  -  C>  x  =  *R  +  +  Xs>  then  =  F(xu+XS)^  Where  V  Xs 

satisfy  (2.15). 


Proof!  This  is  obvious  from  Theorem  1  and  Lemma  2.2. 

The  conditions  of  Theorem  2.2  imply  that  equations  (2.15) 
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(the  determining  equations)  are  equivalent  to  the  equation  Ax  =  Nx 
for  the  solutions  x  =  x^  +  x^  +  x^  with  |x^|  <  d,  |x^+Xg|  <  c. 

It  is  of  interest  to  determine  for  given  c,  d  conditions  on  A 
which  will  ensure  that  the  range  of  the  projection  operator  S  is 
always  finite  dimensional.  Such  a  result  is  contained  in'. 


Theorem  2 .3 .  Suppose  X  is  a  Hilbert  space,  A  satisfies  (Hi)  and 
the  set  V  =  {x  e  (i-U)^(A)!  |Ax|  <  1}  is  relatively  compact. 
Then,  for  any  k  >  0,  there  is  a  projection  S!  X  -»  X  with  finite 
dimensional  range  and  satisfying  (H2),  (H5).  Consequently,  if  (h4) 
is  satisfied,  then,  for  any  positive  constants  c,  d  there  always 
exists  an  S  such  that  (2.l4)  is  satisfied  and,  therefore,  the  con¬ 
clusions  of  Theorem  2.2  hold. 


Proof I  If  we  call  B,  the  unit  ball  in  Z_,  we  see  that  M(B.)  =  V. 
-  1  E*  d 

and  then  M  is  a  compact  (and  hence  bounded)  operator.  Since  V 
is  totally  bounded,  for  any  k  >  0,  there  is  a  positive  integer  n(k) 
and  a  sequence  {x.}  C  V,  j  =  l,2,...,n(k)  such  that  for  any  z  in 
B,,  there  is  an  x.  with  |Mz-x.|  <  k.  Let  S  =  S(k)!  X  -*  X  be  the 
projection  whose  range  coincides  with  the  finite  dimensional  subspace 
spanned  by  {x.. )  and  whose  null  space  is  [x^,  . .  ]*^  ©  X^, 

where  the  orthogonal  complement  is  taken  on  X^  Since  an  ortho¬ 
gonal  projection  has  norm  1,  |l-S|  <  a,  where  a  depends  only  on  U. 
Then,  for  each  z  in  B^  we  have:  |  (i-S)Mzj  =  |  (I-S)  (Mz-x.)  |  <  kcc. 

Clearly  US  =  SU  =  0.  Since  UM  =  0,  Lemma  2.3  implies  (H3)  is 
satisfied.  The  remainder  of  the  proof  is  obvious. 


Theorem  2.4.  Suppose  the  conditions  of  Theorem  2.3  are  satisfied  and 


there  is  a  sequence  of  finite  dimensional  orthogonal  projection 

operates  S  !  X_  TT  -*  XT  TT  such  that  X„  CX„  ,  m  =  1.2....  and 
m  X-U  I-U  Sffi  Sffi+1> 

for  any  finite  dimensional  projection  operator  T!  X^  U  ^1  U  an(^ 

any  e  >  0,  there  is  an  m q(c)  >  0  such  that  | (l-Sm)T|  <  8  for 

m  >  m  (g).  For  any  constants  0  <  c  <  d,  there  is  an  integer  m.  >  0 

such  that  each  S  for  m  >  m.  satisfies  (2.l4)  and,  therefore, 
m  —  j  \  j  >  > 

the  conclusion  of  Theorem  2.2  hold. 

ProofI  For  any  k..  satisfying  (2.l4)  there  always  exists  an  S  = 

S(k-^)  by  Theorem  2.3  such  that  (H2),  (H3)  are  satisfied.  Choose 
m^  so  large  that  I  (l-S  )SM|  <  T)  where  |e]  (rj+k^)  < 
min  {[a(c+d)]_1,  d[p(c+d)]”'1'}  for  m>m..  Then 

|(l-Sm)M|  <  |(l-Sm)SM|  +  1  (I-Sj  (I-S)M]  <  Tj  +  1  (l-S)Ml  <  k.  +  tj. 

Therefore,  Sm  satisfies  (H2),  (H3)  with  k  =  k^  +  t).  On 

the  other  hand,  the  choice  of  tj  implies  that  (2.l4)  is  satisfied 

with  k  =  k.  +  tj.  The  remainder  of  the  proof  is  clear. 

Cl 

In  the  case  where  the  determining  equations  are  finite 
dimensional  and  are  obtained  by  the  application  of  the  contraction 
principle  on  the  operator  A  in  Theorem  2.1,  Williams  [1]  has  shown 
that  it  is  possible  to  establish  a  theoretical  connection  with  the  Leray- 
Schauder  theory  and  topological  degree. 
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3.  Methods  of  Solution 

Under  appropriate  hypotheses  on  A,  we  have  seen  in  the  previous 
section  that  the  solution  of  Ax  =  Nx  can  be  reduced  to  the  solution  of 
equations  (2.2)  or  (2.8).  Our  goal  in  this  section  is  to  discuss  in  a  very 
general  manner  various  methods  for  obtaining  approximate  solutions  of  these 
equations.  An  abstract  version  of  this  problem  is  as  follows;  Suppose 
U,  V,  W  are  Banach  spaces  and  B:  D  X  C;  U  X  V  -»  W  are  continuous 

operators.  The  solution  of  equations  (2.2)  or  (2.8)  is  a  special  case  of 
the  following  problem:  find  u  e  U,  v  e  V  such  that  the  system  of  equa¬ 
tions 


(3.1) 


u  -  B(u,v)  =  0 
C(u,v)  =  0 


is  satisfied. 

_  / 

If  B  and  C  have  continuous  Frechet  derivatives,  the  classical 

iterative  method  for  solving  (3-1)  is  Newton’s  Method.  Let  Bu(u,v), 
Cu(u,v),  B^(u,v),  C^(u,v)  denote  the  Frechet  derivatives  of  B(u,v)  C(u,v) 
with  respect  to  u,v,  respectively.  If  the  linear  mapping  u,v) ; 

U  X  V  -»  V  X  W  given  by 

[.#(u,v)](<p,t)  =  [(I-Bu)cp  -  B n|r,  Cucp  +  Cvil 

has  a  bounded  inverse  for  every  (u,v)  in  some  bounded  set,  then  Newton's 
method  is  given  by 


(3.2) 


(VV  given  J 


-lr 


(Vl’Vll  *  <Vvn>  -  tJ(VTn>]  l".-!(V,u>>C(".'TnH  t11  *  °>1’2> 


Conditions  for  the  convergence  of  this  method  as  well  as  error  bounds  may 
be  found  in  Antosiewicz  [3], 

There  are  many  variants  of  Newton’ s  method  that  can  be  used  to  ad¬ 
vantage.  For  example,  we  could  always  evaluate  the  inverse  of  &  at 
(uq,Vq)  in  order  to  eliminate  some  computations.  Another  scheme  we  might 
choose  is 


(Vv0>  given5 


(3.3)  u  .  =  u  -  [I-B  (u  ,v  )]-1[u  -B(u  ,v  )], 

v  '  n+1  n  L  uv  n’  n'  J  L  n  v  n’  n'  ’ 


v  .  =  v  -  [C  (u  ,v  )]"  C(u  ,v  )  (n  =  0,1,2,...). 
n+1  n  L  vN  n’  n'  v  n*  n7  v  9  7  7  Jy 


provided,  of  course,  that  the  inverse  operators  exist.  Alternatively,  we 
could  attempt  the  iteration 


(3.^) 


(uQ,v0)  given, 


Vl  *  B<VTn>i 


.-1 


vn+l  =  Vn  "  C^VVn)  <n  =  °U,2,...). 


Convergence  criteria  and  error  bounds  are  easily  obtained  for  each  of  the 


above  methods. 
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A  case  of  (3»l)  that  occurs  quite  frequently  in  the  applications  is 
when  N  =  £N,  where  6  is  a  small  real  parameter  satisfying  0  <  |e|  <  £q. 

If  H^)  and  H^)  are  satisfied  for  S  =  0  and  if  N  is  locally  Lipschitzian, 
then  we  can  see  that  the  conditions  of  Theorems  2.1  and  2.2  are  satisfied 
for  8^  sufficiently  small.  Therefore,  there  exists  a  function 

fs  v(c)  x  [-e0,eQ] 

such  that  F(xy,e)  is  continuous  in  Xy,£,  and  satisfies  F(x^,0)  =  0. 

Thus  system  (2.1)  has  a  solution  if  there  is  an  >0,  <  E^,  for  which 

the  equation 

(5.5)  (i-E)N(xu  +  F(xu,e))  =  0 

has  a  solution  x-y  =  Xy(£)  for  0  <  |s|  <  £^. 

If  N  has  a  continuous  Frechet  derivative  IT  ,  then  F(xy,£)  has 
continuous  derivatives  in  Xy^e.  Therefore,  if  there  exists  an  x^  such 
that 

a)  (I-E)N(xg)  =  0J 

(5-6) 

b)  (I-E)N  (x^)  has  a  bounded  inverse, 

then  the  implicit  function  theorem  implies  the  existence  of  a  solution  x-^£) 
of  (5*5)  for  |  £ |  <  £^,  £^  <  £q  sufficiently  small,  satisfying  x^(0)  =  x^. 
Thus,  equation  (2.1)  has  a  solution 


x  =  x  (e)  +  f(x  (e),e)  for  |e|  <e 


-  ~r 


A  rather  natural  way  to  calculate  the  solution  x  for  this  special 
case  is  to  use  a  variant  of  namely, 


xn 


VV  un  e  Xl-u'  vn  €  *0  (n  = 


(3-7) 


uo  =  °>  vo  =  V 


un+l  =eMEN(un+vn); 


vn+i  =  V^CCl-EjNjx^l^d^NCv^u^)  (n  =  0,1,2,...) 


xv  U' 


n  n+1' 


It  is  fairly  easy  to  show  that  this  method  converges  for  |c|  sufficiently 
small  ( see  Lazer  [  1]  for  a  special  case) . 

If  N  is  analytic  in  some  region,  then  the  function  F(x^,e)  will 
be  analytic  in  an  appropriate  region,  and  condition  (3«6)  will  imply  the 
solution  x-y(e)  is  analytic  in  some  region.  An  obvious  way  to  obtain 
the  power  series  expansion  of  x  =  x^(S)  +  F(x^(e),6)  in  g  is  to  simply 


=  kS>  ^  ’  Xk  *  VV  “k  £  XI-«!  Vk  6  X’ 


and  try  to  determine  the  coefficients  successively.  These  coefficients 
must  satisfy 


no  -  °>  vo  '  *v 


Vi  ’  VV-'V  V"'VS 


(5-8) 
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( I-E)Nx(  Xy )  vk+2_  ~  U0 9  *  *  *  9  uk*  V0*****Uk^  ^  =  ®  9  •  •)  > 


where  the  f  9  g  are  some  given  polynomial  functions  of  their  arguments. 
k  k 

Condition  (3.6)  implies  that  these  equations  can  he  solved.  This  is  the 

/ 

natural  generalization  of  the  Poincare  procedure  for  obtaining  periodic 
solutions  of  ordinary  differential  equations. 

In  general,  the  method  (5*7)  is  probably  a  more  accurate  way  to 
obtain  the  solution,  but  (3«8)  is  certainly  very  simple. 

We  conclude  this  section  with  a  brief  description  of  the  notion  of 
quasilinearization.  We  have  seen  that  the  equation  Ax  =  Rx  is  equivalent 
to  the  system 


xl-u  '  ♦  xi_„)s 


(I-E)M(xu«I_u)  »  0 


If  we  let  x  =  v,  Xj  =  u,  the  system  becomes 


u-MEN(u+v)  =  0$ 

( I-E)R( u+v)  =  0. 


For  solution,  we  use  the  modified  Rewton  method: 


un+l  =  Un  '  [I-MER^u+vJ]'  [ u„-MEN( u^+v J ] ; 


x'  n  n 


vn+l  =  vn  -  [(l-E)NxCVvn)]-  (l-E)R(un+vn)  J 
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u  ^  =  ME{N  (u  +v  )[u  n-u  ]  +  N(u  +v  )) ! 
n+1  xv  n  n'L  n+1  nJ  v  n  n'  1  * 

(l-E){[N  (u  +v  )][v  , -v  ]  +  N(u  +v  )}  =  0. 

s  '  L  yr  n  n'  L  n+1  n  v  n  n' 

This  last  system  is  equivalent  to 

[A-M  (u+v)]u  ,  =  E[N(u  +v  )-N  (u  +v  )u  It 

L  xv  n  n' J  n+1  L  '  n  n'  xv  n  n'  nJ  ’ 

(I-E)N  (u  +v  )v  =  (l-E)[N(u  +v  )  -  N  (u  +v  )v  ]. 
v  '  xv  n  n'  n+1  v  ' L  '  n  n'  xv  n  n'  nJ 


Thus  whenever  I  =  E,  the  pair  of  equations  above  reduces  to  a  single  equa¬ 
tion.  In  such  an  event,  we  have  quasilinearization  (see  Antosiewicz  [3]). 
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4.  Periodic  Solutions.  To  illustrate  how  to  formulate  specific 
problems  in  the  above  manner,  consider  the  equation 

x  =  Bx  +  ef (t,x)  (4.1) 

where  f(t+2Tr,x)  =  f(t,x)  is  an  n-vector  continuous  in  (t,x),  B 

is  an  n  x  n  constant  matrix  and  S  is  a  small  real  paramter. 

Our  problem  is  to  find  27T-periodic  solutions  of  (4.1). 

To  formulate  this  in  abstract  terms,  let  X  =  Z  be  the 

Banach  space  of  continuous  27T-periodic  functions  with  |x|  = 

sup  |x(t)|  for  x  e  X.  For  any  x  e  X  with  x  continuous  let 
0<t<27T 

(Ax)(t)  =  x(t)  -  Bx(t),  -oo  <  t  <  oo,  (4.2) 


and  for  any  x  e  X,  let 

(Nx)(t)  =  gf(t,x(t)),  -oo  <  t  <  oo.  (4.3) 

Finding  27T-periodic  solutions  of  (4.1)  is  now  equivalent  to 
finding  an  x  e  X  such  that 


(4.4) 


with  A,N  defined  in  (4.2),  (4.3). 

The  set  ^(A)  is  the  set  of  all  27T-periodic  solutions  of  the 
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unperturbed  equation I 


x  -  Bx  =  0. 


(**•5) 


Therefore,  it  is  a  finite  dimensional  subspace  of  X  and  admits  pro¬ 
jection  by  an  operator  U,  tyt(A)  =  X^.  If  J  is  the  set  of  integers 
j  such  that  (B-jl)  is  singular,  e^, ...,e^  are  a  basis  for  the 


eigenspace  of  j  and  x  ~  Z  x^  exp(ikt),  then  Ux  =  Z  x^exp(ikt) 

k=  -oo  ke  J 


U 


where  x^  is  the  projection  of  x^  onto  the  span  of  the  vectors 


e-^,  •  •  • ,  s 


Pb 


The  classical  Fredholm  alternative  states  that  the  equation 


x  -  Bx  =  y,  yeX 


(4.6) 


has  a  solution  x  e  X  if  and  only  if 


27 r 

/  z(t)y(t)dt  =  0 
0 


(M) 


for  all  row-vectors  z  for  which  z'  e  X  ( ’  =  transpose)  and 

z  +  zB=  0.  (4.8) 

Since  the  linear  subspace  spanned  by  such  z'  is  finite 
dimensional,  if  follows  that  &£{k)  =  X^,  for  some  projection  E. 
To  check  (H3),  we  see  that  for  any  x  e  @(A) ,  x  ~  Z  x,  e^^  f 
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ikxkeik"t,  and  hence 


2  k‘ 


1^1 


<  00. 


Then 


|  (l-U)x|  =|  Z  (ikI-B)"1(ikI-B)^elkt|  < 

k^J 

(  Z  I  (ikl-B )"1|2)1/2(Z  I  (ikl-B)x^l2)1/2  <  /  lAxCt)]^)1/2  < 

k^J  0 

<  M-jAx]  for  some  constant  M^. 

Consequently,  the  previous  theory  implies  (if  f(t,x)  has  a  continuous 
first  derivative  with  respect  to  x)  for  any  c  <  d,  there  is  an 
eQ  >  0  such  that  for  any  fixed  u  e  Xy,  ]u]  <  c,  there  is  a  unique 
solution  F(u,e)  €  X^  ^  of 


Av  =  EN(u+v) 


for  0  <  |g|  <  and  F(u,e)  is  continuously  differentiable  in  u,g. 
The  bifurcation  equations  are  then  given  by 


(i-E)N(u  +  F(u,g))  =  0. 


(M) 


Since  'Jl(A)  =  is  finite  dimensional,  it  is  more  con¬ 
venient  to  write  the  bifurcation  equations  in  terms  of  a  variable  in 
some  Euclidean  space.  If  $  =  (cp^,  .  . . ,9^)  is  a  basis  for  91(A) 
then  u  =  $a  for  some  d-vector  a,  and  the  bifurcation  equations 


become 
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G(a,g)  =  (l-E)N(®a+F(®a,e))  =  0. 


(4.io) 


If  f(t,x)  is  analytic  in  x  in  some  region,  then  F(u,g) 
is  analytic  in  u,g .  One  can  show  that  if  there  is  an  a^  such  that 


G(a0,0)  =  0 

(aQ>®)]  /  0 


(4.11) 


then  there  is  a  solution  a(g)  of  (4.10)  for  g  small,  which  is 
analytic  in  g.  One  could,  therefore,  obtain  the  solution  in  powers 
of  g  by  substituting  in  the  original  equation 


x 


00 


Z  X  g 

k=0 


k 

y 


-  »ak  +  V  Tk  €  xi-u 


(4.12) 


and  determining 


so  that  the  equation  (4.1)  is  satisfied. 

/ 


This 


is  the  classical  Poincare  perturbation  procedure  for  obtaining  periodic 
solutions  of  (4.1). 

Of  course,  any  of  the  other  methods  mentioned  before  give 
other  iterative  schemes.  In  particular,  the  analogue  of  (3-7)  f°r 
this  special  case  and  hypothesis  (4.11)  was  given  by  Lazer  [1). 

In  many  applications,  condition  (4.11)  may  not  be  satisfied 
and  therefore  the  above  iterative  procedures  fail.  However,  it  is 
sometimes  possible  to  discuss  the  qualitative  properties  of  the 
determining  equations  (4.10)  without  any  successive  approximations. 
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We  mention  only  the  classical  problem  of  Liapunov  on  the  existence 
of  two  parameter  families  of  periodic  solutions  of  Hamiltonian  sys¬ 
tems.  One  can  show  without  successive  approximations  that  the 
bifurcation  equations  are  identically  zero  (see  Hale  [4]). 

One  point  about  the  above  procedure  that  has  been  neglected 
is  the  precise  determination  of  the  projection  operators  U  and  E. 
There  are  many  ways  to  choose  these  and  each  gives  a  different  form 
for  the  approximating  equations.  We  mention  only  one  way. 

If  0  is  the  n  x  d  matrix  whose  columns  are  a  basis  for 
the  27r-per iodic  solutions  of  (4.1)  and  ¥  is  a  d  x  n  matrix  whose 
rows  are  a  basis  for  the  2jr-periodic  solutions  of  (4.8),  define  the 
nonsingular  d  X  d  matrices  C,D  by 

2ir  2rr 

C  =  /  ®' (t)fl(t)dt,  D  =  /  T(t)Y' (t)dt 

0  0* 


and  projection  operators  U,Q  on  X  by! 


Ux  =  $a 


Qx  =  Fb 


a 


2rr 

/  3>'  (t)x(t)dt 

0 


1  ar 

b  =  D  /  ¥(t)x(t)dt,* 

0 


then 


Xy  =  91(A),  XX_Q  =  ^(A),  E  =  I  -  Q. 


-26- 


To  illustrate  the  application  of  the  above  ideas  to  a  problem 
involving  no  small  parameter,  we  consider  the  following  problem  of 

Lazer  and  Sanchez  [1],  Suppose  x  is  an  n-vector,  X  is  the  space 

2  n 

of  L  n-vector  functions  of  period  2/r,  p  e  X,  H!  E  -»  R  is  a  func¬ 
tion  with  continuous  second  derivatives  satisfying! 

N2I  <  MjjI  <  (|^-)  <  Mjj+1I  <  (N+1)2I  for  all  a  e  E*. 

i  d 

A  matrix  is  greater  (strictly  greater)  than  zero  if  it  is  positive 
(positive  definite).  The  problem  is  to  discuss  the  existence  of  a 
2/r-periodic  solution  of 

x  +  grad  H(x)  =  p(t)  (4.15) 

The  case  n  =  1  was  first  treated  by  Loud  [1]  and  generalized  by 
Lazer  and  Leach  [1],  For  the  n-dimensional  case,  we  have  the  fol¬ 
lowing 

Theorem.  Under  the  above  conditions,  the  equation  (4.15)  has  a  2rr- 
periodic  solution. 

Proof!  We  only  indicate  the  essential  ideas  of  the  proof  and  the 

reader  can  consult  the  original  paper  for  details.  Let  A!  C 

X  -»X  be  defined  by  (Ax)  (t)  =  ,  t  e  (-«,»)  and  N!  X  -» X 

dt 

be  defined  by  Nx(t)  =  -grad  H(x(t))  +  p(t),  t  e  (-oo,«).  Then 
91  (A)  =  Xy  is  the  set  of  constant  functions  in  X. 

The  operator  A  obviously  has  a  bounded  right  inverse  M 
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and  without  loss  in  generality,  we  can  assume  UM  =  0.  Therefore, 
finding  2ir-periodic  solutions  of  (4.15)  is  equivalent  to  Ax  =  Nx 
which  in  turn  is  equivalent  to 


a)  x  =  Ux  +  M(l-U)Nx 

b)  UNx  =  0. 


(4.16) 


If  x 


y  +  z. 


y  e  z  e  Xy,  then  (4.l6)  is 

a)  y  =  M(l-U)N(y+z) 

b)  UN(y+z)  =  0. 


equivalent  to 


(4.!7) 


The  first  step  of  the  proof  is  to  show  that  for  any  z  e  X^, 
there  is  a  unique  solution  F(z)  of  (4.17a).  This  is  accomplished 
in  the  following  manner!  For  any  p  such  that  (i-pM)  has  a  bounded 
inverse,  equation  (4.17a)  is  equivalent  to 


y=  (I-mM)_1M[  (l-U)N(y+z)  -  py]d=f  (l-pM)_:iTy[G(y,z)  . 

Using  methods  very  similar  to  those  of  the  Hilbert- Schmitt 
theory  of  Fredholm  integral  equations  and  for  p  =  ^  (Pjj+q  +  it 

is  then  shown  that  the  operator  (l-pM)-\lG(  • ,  z)  is  a  contraction 
for  each  z.  This  depends  upon  the  fact  that  M  is  completely  con¬ 
tinuous  and  has  a  complete  orthonormal  system  of  eigenvectors  and 
that  the  Gateaux  derivative  G' (x)  of  G  satisfies 
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|G’(x)|  <|  [ (W+l)2  -  N2]. 

In  the  process  of  the  proof  of  the  contraction  one  also  ob¬ 
tains  that  F(z)  is  globally  Lipschitzian  in  z,  with  Lipschitz 
constant  L  and  |F(z)  -  z]  <  La  for  some  constant  a.  Therefore, 
UN(F(z)  +  z)  is  globally  Lipschitz  in  z  with  Lipschitz  constant  K. 

The  next  step  is  to  show  that  the  function  y(z)  =  TJN(F(z)  +  z) 
satisfies  z-y(z)  -»  «  as  jz]  ->  oo,  and  for  any  r  for  which  z«y(z) 

>  0  on  V  =  [Z*.  |Z|  =  r),  there  is  an  e  >  0  such  that  |z  -  ey(z)| 

2  2 

<  |  z  |  =  r  on  V.  This  implies  the  map  z  -  6t(z)  takes  -*  V 
and  a  known  theorem  in  En  implies  there  is  z ^  e  such  that 
z  -  &t(Zq)  =  Zq"  that  is,  t(zq)  =  0.  This  proves  the  theorem. 

We  now  indicate  the  method  of  Cesari  [4]  for  the  justifica¬ 
tion  of  the  Galerkin  procedure  for  finding  periodic  solutions  of 
ordinary  differential  equations.  This  method  applies  to  general  non¬ 
linear  problems  and  involves  choosing  a  projection  operator  S  as 
in  Lemma  2.2. 


Let  X  be  the  Banach  space  of  continuous  2?r-per iodic  func¬ 
tion  with  the  topology  of  uniform  convergence  and  let  A\  2$ (A)  C 


dx 


continuous 


X  -»  X  be  defined  by  (Ax)  (t)  =  ,  -»  <  t  <  ».  For  a  given 

function  g‘.  ExEn->  En,  g(t,x)  locally  Lipschitzian  in  x, 
g(t+27r,x)  =  g(t,x)  for  all  t,x,  let  HI  X  -» X  be  defined  by  (Nx)  (t) 

g  (t ,  X  (t )  )  ,  -oo  <  t  <  oo. 

The  equation 


lift)  =  g(t,n(t)) 


(4.18) 
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has  a  2Tr-periodic  solution  if  and  only  if  there  is  an  x  e  X  such 
that 


Ax  =  Nx. 


(4.19) 


It  is  clear  that  9i(A)  =  X..  is  the  set  of  constant  functions 

1  27r 

in  X  and  U  may  he  chosen  as  (Ux)  (t)  =^-J  x(s)ds,  _oo  <  t  <  oo. 

Also,  <^(A)  =  A  bounded  right  inverse  M  for  A  is  defined 

for  any  z  e  X^y  as  the  unique  solution  of  mean  value  zero  of  the 
equation  (Ax)  (t)  =  x(t)  =  z(t),  -«  <  t  <  «.  Clearly  Ml  Xj  ^  -*  X 
DM  =  0. 

Any  x  e  X  has  a  Fourier  series 


x(t) 


00 

E 

k=-o 


a.  e 
k 


ikt 


ak  ^7 r 


2ir 

f  e-1  x(t)dt. 
0 


For  any  given  integer  m  >  0,  let  S  :  X  X  be  defined  by 


(Smx)(t)  =  E  a  e 
r-  0<]k|<m  k 


ikt 


(4.20) 


One  easily  sees  that  SmM  =  MSm  for  all  m  >  0  and  Lemma  2.2 
implies  that  equation  (4.19)  is  equivalent  to 
y  =  M(l-Pm)N(y+z) 

?m[A(y+z)-N(y+z)]  =  0 

^m  =  U  +  Sm>  y  e  X  ~  >  z  e  x~  • 


I-P 


m 


m 


(4.21) 
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Using  the  definition  of  P  ,  M,  Parseval' s  inequality  and  Schwarz's 
inequality  one  easily  shows  that 


]M(I-P  )x|  <  r(m)]x],  r(m)  =  (2tt)1//2(  Z  k-2)1/2. 

|k|>m 


Since  ir(m)  -*  0  as  m  ->  »,  one  can  use  the  contraction 
principle  to  prove  for  any  c  >  0,  d  >  0  that  the  conditions  of 
Theorem  2.1  are  satisfied  for  m  sufficiently  large.  Thus  for  m 
sufficiently  large,  the  equation  (4.18)  has  a  27r-periodic  solution 

if  and  only  if  there  exists  a  z  e  X  such  that 

P 

m 

Pm(A-N)(z+F(z))  =  0.  (4.22) 


In  words,  the  above  remarks  imply  the  following!  there  is 
always  an  integer  m  such  that  one  can  fix 


z(t) 


z 

|k|<3n 


z,  e 
k 


ikt 


and  determine  a  function  x  =  F(z)  +  z  in  such  a  way  that  the  Fourier 
series  of  the  function  v(t)  =  x(t)  -  g(t,x(t))  contains  only  the 
harmonics  e1^,  |k|  <  m.  The  determining  equations  (4.22)  involve 
the  determination  of  the  z^,  ]k|  <  m,  in  such  a  way  that  the  remaining 
Fourier  coefficients  of  v(t)  vanish. 

Now  to  determine  the-  existence  of  a  periodic  solution  using 
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this  method,  one  must  determine  first  how  large  m  must  he  in  order 
to  apply  the  contraction  principle  and  then  try  to  show  the  existence 
of  a  solution  of  the  determining  equations  (4.22)  even  though  they 
involve  the  function  Fz  which  can  only  he  known  approximately.  One 
usually  proceeds  in  the  following  manner!  first  obtain  an  a  priori 
hound  on  F(z)  as  |F(z)|  <  5  for  some  constant  8^.  Then  try  to 
show  the  equation 

Pm(A-N) (y+z)  =  0 

has  a  solution  for  every  function  y  e  X  ^  with  ]y|  <  5  . 

I-P  m 

m 

To  show  this  latter  property  is  satisfied  one  naturally  looks 
at  the  problem  for  y  =  0, 

Pm(A-N) (z)  =  0, 
th 

which  is  the  m  Galerkin  approximation  to  the  solution.  An  index 
argument  can  sometimes  be  used  to  complete  the  investigation.  For 
some  applications  see  Cesari  [4],  Locker  [1,  2], 

In  the  limited  space  we  have  here,  we  give  a  simpler  applica¬ 
tion  due  to  Knobloch  [1].  Consider  the  system  of  second  order  differ¬ 
ential  equations 

y  +  f(t,y,y)  =  0 

f(t+27r,y,z)  =  f(t,y,z) 

y  =  (y-L,  t  =  (flf . . . ,tm) 


(4.23) 
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and  assume  that  f(t,y, z)  is  hounded  in  (t,y, z)  and  globally 
Lipschitzian  in  y,z  and  there  exist  <  fL  such  that 


f..(t,y,z)  <0  if  yi  <  a 

f^Ct^z)  >  0  if  yi  >  p±,  i=l,2,...m. 


We  assert  that  the  above  hypotheses  imply  that  (4.2^)  has 
at  least  one  27T-periodic  solution. 

In  fact,  repeating  the  same  process  as  above  with  (Ay) (t)  = 
y(t),  (Ny) (t)  =  -f (t,y(t) ,y(t) ),  one  concludes  that  there  is  an  integer 


n^  and  a  function  FI  X  -> X 

0  p  i-p 

no  Qo 


such  that  F  is  continuous  and 


equation  (4.2J)  has  a  27T-periodic  solution  if  and  only  if  there  exists 

a  u  e  X  satisfying  P  (A-N)u  =  0-  These  equations  are  equivalent 
P  no 

no 


gj(a)  d=f  -J2aj  -  ^  /  f(t,u(t)  +  Fu(t),  u(t)  +  F(u(t))dt  =  0 

2tt  0  <  I  <3 1  <  n0 

S0(a)  d=f  f  f(t,u(t)  +  Fu(t) ,  u(t)  +  Fu(t))dt  =  0 


u(t)  =  Z  a  1 

l*l2>0 


ikt 

a,  e  and 


a  (a_n^>  * ' ' f a_i>  &0’ ai> ' • ' ) an^) ' 
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In  the  space  of  vectors  a,  let  us  consider  now  a  box 

around  the  origin  with  faces  parallel  to  the  planes  a.  =  0, 

d 

0  <  |  j|  <  nQ  and  such  that  gj(a),  0  <  |  j|  <  m,  takes  opposite  signs 
on  opposite  faces  to  the  plane  a^  =  0.  For  these  faces  fixed,  move 
the  other  faces  so  far  apart  that  g^a)  takes  opposite  signs  on 
these  latter  faces.  Mirando' s  version  of  the  Brouwer  fixed  point 


theorem  now  implies  there  is  a  zero  of  all  the  functions 
0  <  |J|  <  n  ,  in  the  box.  This  proves  the  result. 


Further  interesting  applications  of  these  techniques  to 


periodic  solutions  of  ordinary  differential  equations  may  be  found 
in  the  papers  of  Mawhin  [1-5]. 
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5.  Eigenvalues  of  Nonlinear  Operators 

Let  X  be  a  Banach  space  and  C:  ^(C)  CX->X  be  a  linear 
operator;  let  also  B:  X  -» X  be  an  operator  such  that  B(0)  =  0  and  B 
has  a  continuous  Frechet  derivative  B* .  Then,  Bx  =  Dx  +  fix,  with  D  = 
B* (0),  fi(0)  =  0,  ft’ (o),  and  in  fact, 

|fix-fiy|  <p(p)|x-y)  for  |  x|  ,  |  y|  <  p, 


where  p(0)  =  0. 

The  eigenvalue  problem  is  the  following;  Does  there  exist  (jl  f  0 
and  x  f  0  such  that  Cx  =  pBx? 

We  say  that  p^  is  a  bifurcation  point  if  for  any  S  >  0  and 
any  p  >  0  there  is  a  p  in  |  P-Pq|  <  S  and  x  /  0,  |  x|  <  p  such  that 
Cx  =  pBx. 

If  we  write  the  equation  as  (C-PqD)x  =  (p-p^)Dx  +  fix  and 
(C  -  PqD)  has  a  bounded  inverse,  then  there  can  be  no  bifurcation  at  p^. 

In  fact,  x  =  (C-p0D)"1[(p-p0)Dx  +  fix]  and  the  contraction  principle  implies 
a  unique  solution  of  this  equation  for  j P-Pq|  <  6,  |x|  <  p  if  S  and 
p  are  small.  Since  x  =  0  is  a  solution,  p^  cannot  be  a  bifurcation 
point.  Therefore,  to  have  a  bifurcation,  (C-p^D)  must  not  have  a 
bounded  inverse. 

If  we  call  A  =  C-p^D,  N  =  (p-p^)D  +  pfi  we  are  going  to  assume 
as  before  that  $ft(A)  =  X^,  i^(A)  =  Xg  for  some  projections  U  and  E, 
and  that  A  has  a  bounded  right  inverse  M,  with  UM  =  0.  For  example, 
if  C  =  I  and  B  is  completely  continuous,  this  hypothesis  is  always 
satisfied.  For  this  latter  situation  (C  =  I  and  B  completely  continuous) 
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Krasnoselskii  [1]  has  proved  by  using  rotation  of  a  vector  field  that  there 
exists  always  a  bifurcation  at  if  the  generalized  eigenspace  of  A 

is  of  odd  dimension.  We  show  that  the  techniques  mentioned  before  can  be 
used  in  the  more  general  situation  to  attack  this  problem,  but  will  only 
be  able  to  conclude  that  bifurcation  exists  (under  an  additional  weak 
hypothesis  on  A)  if  dim  91(A)  =  1. 

As  before,  the  equation  Cx  =  pBx  is  equivalent  to  Ax  =  Nx,  and 
hence  it  is  equivalent  to  the  system 

(5-1)  a)  x  =Ux  +  MERx 

b)  (i-E)Nx  =  0. 

If  we  call  x  =  y+z,  y  e  Xj_^,  z  e  ( 5 - 1)  can  be  written  as 

(5.2)  a)  y  =  MEN(y+z) 

b)  (l-E)E(y+z)  =  0. 

For  any  positive  constants  c  and  d  let 

V(c)  =  {x  e  XyS  |  x|  <  c),  j^(d)  =  {x  e  X!  |  x|  <  d,  Ux  =  0). 

Theorem  5.1.  There  exist  eQ  >  0,  dQ  >  0  such  that  for  every  0  <  c, 

d  5  d0>  there  exists  a  unique  continuous  function  Ft  V(c)  x  {|p-po|  < 

eQ)  -»  X  such  that,  F(0,p)  =  0, 

F(x,p)  €  ^(d)  for  x  e  V(c)  and  |  p-p  |  <  e0 


and 
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F(x,p)  =  MER(F(x,p)+x),  |p-pQ|  <  eQ. 

Furthermore,  x  =  x  +  F(x,  p)  satisfies  Cx  =  pBx  if  and  only  if 
x  satisfies 

( I-E)N(x+F(x,  p)  )  =  0 

The  function  F(x,p)  has  a  continuous  first  derivative  with  res¬ 
pect  to  x,p  and  there  is  a  constant  a  such  that  |F^(0,p)|  <  a|  p-p  |  . 

x  u 

Phooft  From  the  definition  of  R,  we  have  R(o)  =  0  and  |  Rx-Ry|  < 

|n|  P(c+d)|x-y|  +  |  p-pj  |  D|  |  x-y|  <[|p|P(2d0)  +  |  D|  |  P-PQ|  ]|  x-y|  for 
I  x|  <  c,  1  y|  <  d,  |  xj  +  |  y)  <  c+d  <  2dQ,  |  p-pQ|  <  eQ. 

The  contradiction  principle  gives  the  result  that  there  exist 
dQ  >  0,  eQ  >  0  and  F(x,p),  x  e  V(c),  |  p-pQ|  <  eQ  such  that 

F(  x,  p)  =  MER(F(x,p)  +  x). 

The  continuity  and  differentiability  of  F(x,  p)  also  is  obtained 
from  the  contraction  principle  and  the  differentiability  of  R  with  res¬ 
pect  to  x  and  p.  The  fact  that  x  must  satisfy  the  bifurcation  equa¬ 
tion  is_  clear. 

To  prove  the  last  assertion  of  the  theorem  observe  the  equality 


F(x,  p)  =  MER(F(x,  p)  +  x) 

M-)  =  MER  (F(x,p)  +  x)  (FJL(x,p)  +  I) 

X  x  X 


implies  that 
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But  N  =  (p-p  )D  +  ft’  (x),  with  ft' (o)  =  0  and  so 
X  u 

*1(0,1!)  =  MB(n-n  )D(l^(0,ji)  +  I), 

X  X 

and  then  if  |  p-p^J  is  small,  we  have  F* 

*l(0,U)  =  (h-Hq)T(p), 

X 

where  T((a)  is  a  hounded  linear  operator  for  which  we  can  find  a  hound  in¬ 
dependent  of  p  if  | P-Pq|  is  small. 

From  the  above  theorem,  the  problem  remaining  is  to  solve  the  bi¬ 
furcation  equation.  Suppose  now  that  the  subspaces  9J(A),  and  (I-E)D  91(A) 
are  one  dimensional.  For  C  =  I  and  D  completely  continuous,  this  latter 
condition  is  not  a  restriction.  For  C  a  self  adjoint  partial  differential 
operator,  D  =  I  and  X  one  of  the  usual  Sobolev  spaces,  it  is  also  not 
a  restriction. 

Let  cp  and  |  be  a  basis  for  91(A)  and  (I-E)D9J(A),  respectively. 

rv/ 

If  x  =  acp,  then  the  bifurcation  equations  are: 

0  =  ( I-E)N(x+F(x,  p)  )  =  (l-E)N(aCp+F(aCp,p)) 

=  (l-E)[(p-pQ)D(acp+F(acp,p))  +  pft(acp+F(acp,p) )  ] 

=  (t!-H0)PM<  +  (p-P0)(I-E)DF(a£p,p)  +  p(  I-E)ft(  acp+F(  acp,  p) ) , 
where  p\| ;  =  (i-E)Dcp,  p  /  0.  If  we  let 
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(l-E)DF(acp,p)  =  *r(a,p) 
p(l-E)ft(a<P+F(asp,p))  =  ^6(a,p) 

the  bifurcation  equations  are  equivalent  to 

(5.3)  (p-pQ)pa  +  (p-p0)F(a,p)  +  5(a,p)  =  0 

From  Theorem  5.1,  there  are  constants  k,  a,  and  a  continuous  func¬ 
tion  r(r),  r  >  0,  Tr(0 )  =0,  such  that  the  following  inequalities  holds 

|F(a,n)|  <  k|  F(acp,  p)|  <  k[o:|  acp|  |  +  r(|  acp|  )*|  acp|  ], 

|S(a,q)|  <k[p(|acp|  +  |  F(sp,p)|)][|  acp|  +  |  F(  atp,  m-)]] 

<kp(|aq>|  +  |F(acp,|i)|  )[|  acp|  +  ajaip|  •  |  p-pQj  +  1f(|  atp|  )|  acp|  ] 

From  these  properties  we  have  the  bifurcation  equations  are  equivalent  to 
(for  any  a  f  0) 


4  -  un  -  (tt-io  il 

0  0  ap  ap 

where  and  ^  are  continuous  as  a  =  0  and  go  to  zero  as 

a  -»0,  p  -»  pQ  (actually  continuously  differentiable).  The  implicit  func¬ 
tion  theorem  implies  there  exists  a  p(a)  satisfying  the  bifurcation  equa¬ 
tion  for  a  and  |  p-pQ|  small  enough.  Also  p(a)  -»  pq  as  a  -»0.  We 
thus  have  the  following! 

Theorem  5.2,  Under  the  above  hypothesis  on  the  simple  eigenvalue  pQ, 
there  always  is  a  bifurcation  at  p  =  p^. 

We  can  find  the  solution  iteratively  ass 
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yo  =  °>  = 

yn+l  =  +  arp) 

r  (a)  5  fa) 

,  \  n '  *  nv  J 


where 


(i-E)Dcp  =  il/p 

(l-E)Dyn  =  *rn(a) 

Pn(l-E)fl(acp+yn)  =  t5n(a). 

These  ideas  are  closely  related  to  the  paper  of  H.  Keller  [1],  where  the 
author  studies  the  equation* 

Lu  +  Xg(X,x,\ i)  *  0,  x  €  D,  Bu  =  0,  x  e  bn, 

L  elliptic  operator  of  order  2m,  B  linear  of  order  m. 

Another  application  of  these  ideas  can  be  seen  in  the  work  of  J.  B. 
Keller  [1].  To  see  this^  consider  the  equation 

a)  [I(t)ut]t  +  f(u,t,\)  =0  ti  <  t  <  t2 

(5.4)  a^Cti)  +  PjuCtp  =  0 

b) 

a2ut(-t2)  +  p2u(t2)  =0,  ax  f  0,  a2  +  p2  0. 

where  u(t)  is  a  scalar  function  and  l(t)  is  positive  and  continuously 
differentiable  while  f  has  derivatives  with  respect  to  \f  u  up  through 
order  three  which  are  continuous  in  some  rectangle  u^  <  u  <  Ug, 


-4o- 


S  ^  S  ^2}  "^1  —  ^  —  ^2* 

Assume  (5.4)  has  a  solution  uQ(t,X)  for  X^  <  X  <  Xg,  u^  < 
u^(t,X)  <  Ug.  Also,  assume  the  linear  variational  problem  at  X^  has  a 
nontrivial  solution,  i.e.,  there  exists  a  solution  cp  of 

(IVt  *  fu(uo(t’V’t>Vv  =  0 

(5.5)  a1vt(t1)  +  *  0 

a2vt(t2)  +  p2v(t2)  =  0. 

Then  necessarily  cp(t1)  =  0  (since  cp^)  =  0  would  imply  9t(t1)  =  0 
and,  by  uniqueness,  cp(t)  would  be  zero).  Also  the  first  boundary  condi¬ 
tion  says  that  cp^(t-^)  can  be  determined  from  cp(t^),  and  so  there  is  only 

one  linearly  independent  solution  of  (5.5). 

^2  2 

Let  us  normalize  it  by  putting  /  cp  dt  =  1.  The  problem  is  to 

tl 

find  conditions  under  which  there  are  solutions  of  (5.4)  near  X  =  XQ 
other  than  the  given  solution  u^(t,X). 

To  phrase  this  problem  in  the  language  of  bifurcation  of  a  non¬ 
linear  operator,  let  X  be  the  Banach  space  of  the  continuously  dif¬ 
ferentiable  functions  on  [tptg]  satisfying  the  boundary  conditions  (5.4b) 
with  the  topology  of  uniform  convergence  of  the  function  and  its  deriva 
tive.  Let 


C:  0(C)  C  X  -*X,  (Cu)(t)  =  (Iut)t,  tx  <  t  <  t2, 

t,  <  t  <  tg. 


(5.6) 


B^:  x  -> x, 


(Bxu)(t)  =  -f(u(t),t,X), 
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The  above  hypotheses  imply  that  Uq(*,X)  e  X  is  a  solution  of 

(5-7)  Cu  =  B^u 

for  X^  <  X  <  X2  and,  for  ^  =  Xq,  the  equation 

(5.8)  (C-B«^(u0(.^d))u  =  0 

has  a  one  dimensional  subspace  of  solutions  spanned  by  cp.  The  problem  is 
to  find  conditions  on  B..  which  will  ensure  that  (5*7)  has  solution 
/  Uq(*,X)  for  each  X  in  a  neighborhood  of  X  =  X^. 

This  is  not  exactly  a  special  ase  of  the  problem  discussed  at  the 
beginning  of  this  section,  but  the  ideas  used  there  are  easily  adapted  to 
this  situation.  If  we  let  u  =  u^(-,X)  +  v,  then  equation  (5.8)  becomes 

(5.9)  [C-B^(u0(*,X0))]v  =  Bx(uQ(.,X)+v)  -  Bx(uQ(.,X))  -  (uQ(  *,  XQ))v 
or  v  satisfies  the  differential  equation 

(5.10)  (Ivt)t  +  fu(u0(t,X0),t,XQ)v  =  -f(uQ(t,X)+v,t,X))  +  f(uQ(t,X),t,X) 

t  v- 

and  the  boundary  conditions  (5.4b).  If 


Nv  .  Bx(u0(-,X)+v)  -  -  B!  (u0(.,X0))v 

then  (5.9)  is  equivalent  to  Ax  =  Nx.  Furthermore,  9f(A)  »  X,  = 
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[u  e  X:  u  -  acp  for  some  constant 
The  function  Rv  has  the  form 


a);  R(A)  =  XT  TT  =  {u  e  Xj  /  cpu  =  0). 


Nv  =  [f  u^..  +  f  ,](X  -X)v  -  f 
L  uu  OX  v  n  ' 


uX  '  0 


where  the  functions  are  evaluated  at  some  intermediate  points  between 
Uo(‘,X)jX  and  u^(*,Xq),  X^.  Therefore,  N  has  the  same  form  as  the  N 
discussed  at  the  beginning  of  this  section.  The  self-adjointness  of  the  dif 
ferential  operator  implies  Theorem  5*1  may  be  applied  to  obtain  the  bifur¬ 
cation  equations  in  the  form  of  (5.3) .  Calculating  the  constants  in  (5-3) 
for  this  particular  case,  the  bifurcation  equations  are 


a[p(X-X^)  +  Ta  +  •••]  =0 

(5-11)  »  -  *  fux)^dt 

t2  5 

r  =  /  f  cp^dt 

,  UUT 

tl 

where  all  functions  in  the  integrand  are  evaluated  at  u^(t, X^),t, X^. 

Consequently,  if  p  ^  0,  T  ^  0,  the  implicit  function  theorem  im¬ 
plies  there  is  a  solution  of  (5.4)  different  from  u^(*,X)  in  a  neighbor¬ 
hood  of  X  =  Xq. 

An  interesting  paper  on  determining  the  stability  of  bifurcating 
solutions  by  using  Leray-Schander  degree  has  been  written  by  Sattinger  [ 1] . 
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6.  Analytic  Systems,  Regular  and  Irregular  Singular  Points. 

In  this  section,  we  indicate  the  manner  in  which  some  of  the  classi¬ 
cal  results  on  the  existence  of  analytic  solutions  of  linear  differential 
systems  with  or  without  singular  points  can  be  obtained  using  the  above 
procedure.  The  presentation  follows  closely  the  paper  of  Harris,  Sibuya 
and  Weinberg  [1]*  The  exploitation  of  the  procedure  for  more  complicated 
solutions  has  not  been  but  should  be  undertaken. 

Suppose  D  =  diag(d, , . . . ,d  )  where  each  d.  is  a  positive  integer, 
D  d  dn  ° 

t  =(t  ,...,t  ),  B(t)  is  an  nxn  matrix  analytic  at  t  =  0,  with  an 

absolutely  convergent  power  series  for  1 1|  <5. 

Our  problem  is  to  find  necessary  and  sufficient  conditions  for  the 
existence  of  analytic  solutions  of  the  equation 


(6.1) 


tDy  =  B(t)y 


To  phrase  this  problem  in  the  language  of  Section  2,  let  X  be  the 
set  of  all  n-vector  valued  functions  x(t)  whose  components  have  ab¬ 
solutely  convergent  power  series  expansions  for  |  t|  <6,  5  >  0.  For 
00 

v*  k 

x  e  X,  x(t)  =  Z  x  t  ,  define 
k=0  * 

M  =  2  I  xkl  6k* 

k=0 

With  this  norm,  X  is  a  Banach  space. 

For  any  x  e  X,  let  Aj  X  -> X  be  defined  by 
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Then,  9?(A)  =  X^,  (Ux)(t)  =  xQ,  |t|  <6,  is  the  set  of  constant 

functions  in  X  and  ^(A)  =  (y  e  X:  y  =  t  x,  x  in  X}. 

If  we  define  X  -» X  by 


d  -1  v  d  -1  , 

x)(t)  =  col(  Z  x  t  En  x  t  ),  1 1|  <  5, 

k=0  1  k=0  n 


where  x^  =  col(xk  , ...,x^  )  then  ^  is  a  projection  on  X  and  &  (A)  = 

1  n 

XT  Q  .  Furthermore,  a  right  inverse  M  of  A  on  i^(A)  is  given  by 
1  T> 

oo  k+1  oo  +k 

(Mx) (t)  =,  (  2  ^  -jj!  M  d) 

k=0  1  1  k=0  n  n 

l-d1  l-5n 

and  |  Mx]  <  max(  5  ,...,6  )  |  x|  .  Therefore  M  is  bounded. 

If  Ns  X  -»X  is  defined  by  (Nx)(t)  =  B(t)x(t),  1 1|  <5,  then  N 
is  a  bounded  linear  operator  and  finding  a  solution  of  (6.1)  in  X  is  equi¬ 


valent  to  finding  an  x  e  X  such  that 


(6.2) 


Ax  =  Nx 


From  Lemma  2.1,  (6.2)  is  equivalent  to: 
a)  y  =  M(I-QD)N(y+z) 

(6-3) 

*0  QjjN(  y+z)  =0,  ye  XI_U,  z  e  Xy 

If  the  linear  operator  M(l-Q^)N  are  a  contraction  operator,  then 
(6.3a)  would  have  a  unique  solution  y  =  Fz  where  F:  Xy  -*  X-j.  y  is  a 
bounded  linear  operator.  A  solution  of  (6.1)  would  then  be  determined  by 
solving  the  finite  set  of  linear  equations  Q,dN(Fz+z)  =  0.  Unfortunately, 
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this  operator  may  not  be  a  contraction.  To  obtain  a  contraction  we  pro¬ 
ceed  as  in  Lemma  2.3.  For  r  >  0  an  integer,  let  S  :  X  -»  X  be  the  pro- 

r+l  k  r 

jection  operator  defined  by  (S  x)(t)  =  2  x,t  ,  1 1|  <5. 

r  k=l 

If  =  I-U-S^,  then  it  is  easy  to  see  that 


RrM(i-y  = 


and,  furthermore,  from  Lemma  2.3  that  (6.2)  is  equivalent  to 


a)  v  =  M(l-QD+rI)N(v+u) 

(6.4) 

b)  VrI(A"N)(v+u)  =  0,  v  e  ^  ,  u  e  Xj._R  . 

r  r 

It  also  follows  immediately  from  the  definition  of  M  that  there 
is  a  (3  >  0  such  that 


(6.5) 


lM(l"Vrl)l  <P/(r+l) 


Consequently,  for  r  sufficiently  large,  the  operator  M(l-Q^  ) 

is  a  contraction.  Therefore,  for  any  fixed  polynomial  u  of  degree  r, 

there  is  a  unique  solution  Fu  of  (6.4a)  with  FJ  Xj  R  -» X-  a  bounded 

r  r 

linear  operator.  This  implies  that  the  equation  (6.1)  has  an  analytic  solu¬ 
tion  if  and  only  if  the  polynomial  u  can  be  chosen  so  that 


(6.6) 


%+ri(A-N)(Fu+u)  =  0 
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Equation  (6.6)  is  a  polynomial  vector  equation  with  the  jth  com¬ 
ponent  of  degree  d.+r-l  for  the  n(r+l)  coefficients  of  the  vector 

3  n 

polynomial  u  of  degree  rj  that  is,  these  are  at  most  nr  +  Z  d. 

d-i 0 

linear  homogeneous  equations  for  nr+n  unknowns.  Therefore,  there  are 
n 

at  least  n  -  Ed.  linearly  independent  solutions.  This  proves  the  follow- 
ing  theorem  of  Perron-Lettenmeyr: 


n 


Theorem  6.1.  Equation  (6.1)  has  at  least  n  -  Zd.  linearly  independent 

0=1  J 

analytic  solutions  for  )  t|  <6. 

Now  let  us  consider  the  case  of  the  regular  singular  point  (d.  =  1 

cl 

for  all  j).  In  this  case,  QD+rI=  ^r+1jx  =  u+sr5  that  is>  %+Yjx 
represents  the  vector  polynomial  of  degree  r  consisting  of  the  first 
(r+l)  terms  of  the  expansion  of  x.  Theorem  6.1  gives  no  information  in 
this  case,  but  we  still  may  draw  some  conclusions.  The  function  (A-N)Fu 
has  a  power  series  expansion  beginning  with  terms  of  degree  (r+l).  There¬ 
fore,  equation  (6.6)  is  equivalent  to 


(6.7) 


(U+Sr)(A-N)u  =  0 


x  ,  V 

If  u  =  Z  u,t  ,  B(t)  =  Z  B  ' t  ,  |t|  <6,  then  equation  (6.7)  is 
k=0  k=0 

equivalent  to  the  system  of  linear  equations 


(6.8) 


(“-b0K 


k 

Z  B.u,  .,  k  =  0,1,. ..,r,  u  1  =  0. 
j=l  J  -1 


Equations  (6.8)  are  the  equations  for  the  first  r+l  coefficients 
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of  a  formal  solution  u(t)  =  Z  at  of 

k=0  * 


(6.9) 


ttl  s  B(t)u 


Consequently  we  have  proved: 

Theorem  6.2.  Any  formal  solution  of  (6. 9)  is  an  analytic  solution. 

Let  us  now  discuss  the  existence  of  solutions  of  (6.9)  of  the  form 
u(t)  =  t  v(t)  where  v  e  X.  Then 

tv  =  [B(t)  -  XI]v(t) 


Let  n  be  the  dimension  of  and  let  N  be  the  number 

X  ox 

of  linearly  independent  solutions  u  of  the  above  form.  The  determining 


equations  corresponding  to  (6.8)  are 


[(k+X)l-BQ]vk  =  Z  B.vk_.,  k  =  0,1,. ..,r 
0  — 3_ 


if  v(t)  =  E  v  t  ,  Clearly  N  <  n  +  n  +...  and  N  >N  ... 

^ _ Q  X  ~  X  X  A,+l 

We  now  show  that  N  >  n  .  Without  loss  in  generality,  we  can  as- 

sume  X  =  0  and  is  in  Jordan  canonical  form.  Each  vector  in  ^I(Bq) 

corresponds  to  a  zero  row  of  Bq  and  the  corresponding  d  .  may  be  reduced 

iJ 

to  zero.  Therefore,  Theorem  6.1  implies  at  least  n  -  (u-Uq)  =  nQ  solu¬ 
tions  analytic  at  t  =  0.  This  proves 
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Theorem  6.3.  max(  n^,  n^+1, . . . )  <  <  n^  +  n^+1  +  ...  .  In  particular,  if 

Bq  =  0,  that  is,  system  (6.9)  has  no  singular  point,  then  there  is  funda¬ 
mental  system  of  solutions  analytic  at  t  =  0.  If  B^-kl  is  a  nonsingular 
for  each  positive  integer  k,  then  Nq  =  n^  and,  in  fact,  for  any  u^  e 
9KBq)  there  is  a  solution  u(t)  of  (6.9)  with  u(0)  =  u^. 

Similar  results  can  be  obtained  by  the  same  methods  for  nonlinear 
systems  tDx(t)  =  f(t,x(t))  (see  Harris  [1]). 

A  more  complete  discussion  of  (6.9)  is  obtained  by  trying  to  find 

G 

a  fundamental  matrix  solution  of  (6.9)  of  the  form  U(t)  =  P(t)t  ,  where 
G  is  a  constant  n  X  n  matrix  and  P(t)  is  an  n  X  n  matrix  analytic 


at  t  =  0.  We  indicate  how  some  information  is  easily  obtained  from 

B 

Theorem  6.1.  If  G  =  B^,  then  U  =  Pt  ^  a  matrix  solution  of  (6.9)  im¬ 
plies 


(6.10)  tP  =  BqP  -  PBq  +  Bx(t)P,  B^O)  =  0. 

2 

If  we  look  at  this  equation  as  an  equation  for  the  n  components 
p  of  P,  then 

(6.11)  tp  =  C0P  +  C^(t)p,  cL(0)  =  0. 

If  no  two  eigenvalues  of  B^  differ  by  a  positive  integer,  then 
no  eigenvalue  of  is  a  positive  integer  and  Theorem  6.5  implies  there 

are  as  many  linearly  independent  solutions  of  (6.11)  as  dim  ^(CQ). 
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Furthermore,  any  e  «o0)  yields  a  solution  p  of  (6.11)  analytic 
at  t  =  0  and  p(0)  =  p^.  From  the  definition  of  p  and  Cq  the  p 
corresponding  to  P  =  I  is  in  91(Cq).  Therefore,  (6.10)  has  a  solution 
P  analytic  at  t  =  0  and  P(0)  =  I.  This  gives  the  classical 

Theorem  6.4.  If  no  two  eigenvalues  of  Bq  differ  by  a  positive  integer, 

then  there  is  a  fundamental  solution  U  of  (6.9)  of  the  form  U(t)  = 

B 

P(t)t  ^  where  P  is  analytic  at  t  =  0  and  P(0)  =  I. 

Golomb  [1]  has  given  some  interesting  results  on  the  expansion  of 
fundamental  solutions  of  (6.1)  even  in  the  case  of  irregular  singular  points. 
The  approach  applies  to  even  more  general  equations  and  the  ideas  are  very 
closely  related  to  those  of  Section  2. 
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7.  Partial  Differential  Equations.  In  this  section,  we  indicate 
how  the  above  ideas  have  been  used  to  discuss  some  particular  problems 
in  partial  differential  equations.  We  do  not  treat  the  most  general 
situation,  but  try  to  bring  out  the  main  ideas  and  give  appropriate 
references  for  more  sophisticated  results. 

Consider  the  wave  equation  in  one  dimension! 

u,>  -  u  =  eg(t,x,u,u, ,u  ),  where 
tt  xx  ’  >  3  t’  x'>  (71) 

g(tjxju,P,(l)  =  -g(‘t,-x,-u,-p,q)  =  g(t+2ir,x,u,p,q)  =  g(t,x+2ir,u,p,q) 

2 

for  all  (t,x,u,p,q)  e  fl(p)  =  {(t,x)  e  R  =  (-oo,oo)x(-oo,oo);  |u|  + 

IpI  +  kl  <  p} .  The  problem  is  to  determine  conditions  on  g  which 
will  ensure  that  equation  (7*1)  has  a  solution  u(t,x)  satisfying 

u(t+2ir,x)  =  u(t,x+27r)  =  u(t,x)  =  -u(trx)^  (t, x)  e  R2.  (7-2) 

Due  to  the  form  of  g  any  solution  satisfying  (7-2)  will  also  satisfy 

u(t,0)  =  u (t ,tt)  —  0,  t  e  (-<»,<»),  x  e  [0,7 T].  (7-3) 

Therefore,  we  are  solving  for  27T-per iodic  solutions  in  t  of  the 
vibrating  string  with  ends  fixed  at  x  =  0  and  x  =  tt. 

Let  R  =  (-oo, oo)  and  define 
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2 

T  =  C<P -  R  ->  R}  9  continuous,  cp(t+2Tr,x)  =  cp(t,x+2Tr)  =  cp(t,x) 

=  -9(VX)K 


W  =  {9 1  R  -»  R,’  cp  continuous,  cp(t,x)  =  p(x+t)  -  p(-x+t), 

p(t+27r)  =  p(t)}; 


*  2 

C-,  =  {9!  R  ->  Rj  cp  and  cp  continuous},' 

X 


=  (9*  R  -»  R,"  9  continuous  together  with  derivatives  up  through 
order  2},* 


x  =  t  n  c2;  z  =  t  n  c*. 


All  topologies  in  these  spaces  are  those  of  uniform  convergence,  and 


norms  will  be  desiginated  by 


lT>  I  lw>  I  *li>  I  1 2>  l  ip  I  \z> 


respectively.  Define  operators  A,W,  and  Q,  by 


a:  3>(k)  c  x  -» z;  (Au)(t,x)  =  utt  -  u^; 


n:  X-»z;  (Nu)(t,x)  =  eg(t,x,u(t,x),  ut(t,x),  ux(t,x));  (7-4) 

,  2tt  2tt 

Q,:t-^T;  (Qcp)  (t,x)  =  «-  /  9(s,x+t-s)ds  -  «-  /  9(3  -x+t-s)ds. 

"  0  w  0 


It  is  easy  to  see  that  Q  is  a  projection  of  T  onto  T~  =  W.  We 

Q 

can  also  consider  Q  to  be  a  projection  on  X  and  Z,  and  as  such, 
it  is  clear  that  3R(A)  =  X^.  Finding  a  solution  of  (7-1)  is  equivalent 
to  finding  a  u  €  Qj(k)  such  that 
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Au  =  Ku.  (7-5) 


The  following  lemmas  give  more  detailed  properties  of  the 
operators  Q  and  A. 


Lemma  7 For  any  ip  e  I,  the  following  are  equivalent! 

i)  <P  e 

2t T 

ii)  /  cp(s,y-s)ds  =  0  for  all  y; 

0 

r  ir 

iii)  /  /  cp (tJx)Tf(t>x)dxdt  =  0  for  all  X  in  T  =  W. 

0  0  Q 


Proof I 

(i)  — >  (ii)!  If  cp  e  then  (Qfp)  (t,x)  =  0  implies 

(Qfp)(t,-t)  =  0;  i.e. 


27 r  27 T 

J  cp(s,-s)ds  =  /  cp(s,2t-s)ds  for  all  t. 
0  0 


Therefore 


4ir  2t r  2tt  4tt 

0  =  /  (QfP)  (t,  -t)dt  =  2/  cp(s,-s)ds  -jj-/  (/  cp(s,2t-s)dt)ds 

0  0  ^00 

27T  27T  2tt 

=  2/  cp(s,-s)ds  -  -  /  (/  qp(s,2t-s)dt)ds 

0  T  0  0 

2rr  2tT  1  27T-S 

=  2/  cp(s,-s)ds  -  /  (__  /  cp(s,u)du)ds 

0  0  ^  -s 

27T 

=  2  /  cp(s,-s)ds 
0 
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since  cp(s,-u)  =  -cp(s,u).  This  proves  (ii). 

If  (ii)  is  satisfied,  then  clearly  (i)  is  satisfied.  To 
prove  the  equivalence  of  (ii)  and  (iii),  suppose  y(t,x)  = 
p(x+t)  -  p(-x+t)  and  observe  that 


27T  7T  2tT  t+7T 

/  /  <P(t,x)[p(x+t)  -  p(-x+t)]dxdt  =  /  /  <p(t,y-t)p(y)dydt 

0  0  0  t-ir 


7T 


-  / 


-7 T 


2ir 

I  9 (t,y-t)dt 
0 


p(y)dy 


Lemma  7.2.  ^?(A)  =  ZT  A 

— — —  1-lqJ 


m:  z 


I-Q 


and  A  has  a  bounded  right  inverse 


Proof:  For  any  cp  e  Z,  consider  the  function 

j  *Jj  0 

U(t,x)  =  i  /  /  cp(0,|)d£d0, 

0  x-t+0 


which  belongs  to  Cg,  and  satisfies  Utt  -  Uxx  =  Cp,  U(t,x+2ir)  =  U(t,x), 

U(t,  -x)  =  -U(t,x) .  To  prove  that  ^?(A)  =  ZT  ,  it  is  sufficient  to 

I“Q 

show  that  U(t,x)  =  U(t+2-ir,x)  if  and  only  if  cp  e  ZT  A.  A  straight- 

I-Q 

forward  calculation  using  the  fact  that  cp(t,-x)  =  -cp(t,x)  gives 


U(t+27T,X)  =  J  / 


t+27T 


~X+27T+t-0 

I  9(0,l)d| 

X-27T-t+0 


d0 


£  J 
0 


t+27T' 


X+t-0 

/  <P(e,l)d6 

x-t+0 


d0 
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2t  x+t-0  ,  t  x+t-2ir-0 

=  i  /  /  «p(e,Od|d e  +  ifj  cp(e,|)d|de 

0  x-t+0  0  x-t+2ir+0 


2^  x+t-0 

=  i/  /  cp(0,|)d|d0  +  U(t,x). 


0  x-t+0 


Therefore,  u(t+27T,x)  =  U(t,x)  if  and  only  if 


2t t  x+t-0 

♦(t,x)  -  n=  /  /  cp(0,|)d|d0  =  0 

11  0  x-t+0 


for  all  t,x.  This  is  equivalent  to  requiring  that  ^(t^x)  =  0  = 
i|rx(t,x),  for  all  x,t,  and  that  ty  is  zero  at  one  point.  Using  the 
fact  that  cp(t,-x)  =  -cp(t,x),  we  obtain 

it(t,x)  =  (Qfp)(t,x); 

1  27r 

♦x(t,x)  =  (Qfp)(t,x)  +  -/  cp(0,-x+t-0)d0; 
t(t,t)  =  -\|r(t,-t). 

Since  \|/(0,0)  =  0,  \|r(t,x)  =0  is  equivalent,  by  virtue  of  Lemma  7*1, 
to  (Qfp)(t,x)  =  0.  This  proves  that  J$(A)  =  Taking  Mp  = 

(U-QU)cp,  the  remainder  of  the  lemma  is  obvious. 


Lemma  7*2  implies  that  Au  -  Ku  is  equivalent  to 


u  =  Qu  +  M(l-Q)Nu, 
QNu  =  0. 
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We  can  prove  the  following  theorems  by  a  straightforward  application 
of  the  contraction  principle. 

Theorem  7.1.  Suppose  g  has  continuous  first  and  second  derivatives 
with  respect  to  x?u}p^q  in  fi(p)  and  let  0<a<b<p  be  given 
positive  constants.  Then  there  is  an  e-^  >  0  with  the  following 
property!  corresponding  to  each  y  e  =  W  H  Cg,  <  a>  and 

each  e,  ] e|  <  there  is  a  unique  function  r  =  r(r,  e)  in  X, 
or  =  r,  |r|x  <  b,  continuous  together  with  its  first  derivatives  in 
•r,  €,  satisfying  F  (  r,  0 )  =  X,  F^f'^O)  =  I,  and 

rtt  -  rxx  '  eBr  -  e«Kr- 

Theorem  J.2.  Under  the  hypotheses  of  Theorem  7*1,  if  there  is  an 
and  a  function  y(e)  in  a,  |  e|  <  e2>  sucl1 

that 


QMT(r(e),e)  =  0,  (7-6) 

then  r(r(e), e)  is  a  solution  of  (7*1)  for  | e|  <  €p.  Conversely ? 
if  (7-1)  has  a  solution  u(e)  e  X,  continuous  in  e  for  | e|  <  €p, 
|u(e)|x<b,  |Qu(e)|x  <  a,  |  e]  <e2>  then  u(e)  =  r(x(e),e),  where 
y(e)  =  Qu(e)  and  r(y(e),e)  satisfies  (7-6). 

Equations  (7*6)  are  the  bifurcation  equations  and  a  solution 
■)f(e)  of  (7.6)  is  necessary  and  sufficient  (as  described  in  Theorem  7*2) 
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for  the  existence  of  a  solution  of  (7*1)  for  e  small.  A  more 
convenient  form  for  equations  (7.6)  is  obtainable  from  Lemma  7.1- 
More  specifically,  equations  (7-6)  are  equivalent  to 


H(r,  e)  =  0,  r  =  r(r,e);  Where 


2tt  H(r,e)  (y) 


def 


2ir 

5  g(s,y-B,r(s,y-s),r.(s,y-s),r  (s,y-s))ds 

0  t  x 


(7-7) 

(0  <  y  <  2tt) 


The  explicit  formulae  for  H(y,0)  and  H^(y,0)  are 


2  tt 

2?r  h(t,o)  (y)  =  /  g(s,y-s,r(s,y-s),r+(s,y-B),r_(s,y-s))ds;  (7.8) 

0  ^  x 

2tt 


27T[H'(r,0)A](y)  =  / 

'  0 


[guA(s,y.s)  +  gpA(.(s,y-s)  +  gqZ^.(s,y-s)]ds  (7 -9) 


where  A  is  an  arbitrary  element  of  X  and  the  arguments  of  g 
gp, g^  are  the  same  as  in  (7.8). 

The  functions  H,H|.  are  continuous  mappings  from  Xq  into 
the  space  of  continuously  differentiable  27r-periodic  functions. 

Thus  we  can  apply  the  implicit  function  theorem  to  obtain  a  solution 
of  (7-7)  and  thereby  a  solution  of  (7.1).  More  specifically,  if 
there  is  a  y  €  Xrt,  I 1  <  a  such  that 

0  1  0'x 


h(t0,o)  =  0; 

[H^o))-1:  Pl-xQ 


(7.10) 


is  bounded  then  there  will  be  a  solution  of  (7-7)  for  e  small. 
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Remark  7^1 ♦  Since 

(QI£T)(t,x)  =  H(r,e)(x+t)  -  H(y,  e)  (-x+t) , 

we  could  apply  the  implicit  function  theorem  directly  to  ob¬ 

serving  that  the  mean  values  of  the  functions  in  the  domain  of  the 
inverse  of  the  derivative  operator  are  unimportant.  This  remark  is 
convenient  in  the  applications . 

To  see  the  nature  of  the  bifurcation  equations  and  especially 
the  functions  in  (7.8),  (7*9),  we  consider  only  one  example.  Let  us 
examine  the  equation 

utt  -  uxx  =  e|>t+bu+cu5+f(t,x)]  (7.11) 

where  b,c  are  constants,  b  ^  0,  and  f(t,-x)  =  -f(t,x).  The  func¬ 
tions  H(r,0),  H^(r,0),  in  (7.8),  (7*9)  are  given  by 

r( t,x)  =  p(x+t)  -  p(-x+t);  A(t,x)  =  q(x+t)  -  q(-x+t); 

H(t,0) (y)  =  pT(y)  +  [b+3cm(p2)]p(y)  +  cp5 (y)  -  cm(p5)  +  h(y); 
(HJ,(t,0)a)  (y)  =  q' (y)  +  [b+6cm(p2)]q(y)  -  3c[m(p2q)  -  2m(pq)]; 

1  2ir 

h(y)  =  5=  /  f(s,y-s)ds,  0  <  y  <  2tt; 

0  "  “ 

1  27f 

where  1  =  d/dy,  m(cp)  =  — —  /  cp(y)dy  and  we  have  assumed  without  loss 

e~*'  0 
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of  generality  that  m(p)  =  m(q)  =  0.  For  c  sufficiently  small 
one  can  show  there  is  a  p  such  that  the  conditions  (7.I0)  are 
satisfied  and^  thus,  there  is  a  solution  u  e  X  of  (7. 11)  (see  Hale 

[3])- 

Many  more  problems  of  the  above  type  have  been  considered 
for  equation  (7*1)  and  for  even  more  general  equations.  We  have  only 
tried  to  indicate  the  procedure  and  refer  the  interested  reader  to  the 
papers  of  Cesari  [6-9],  Rabinowitz  [1],  Vejvoda  [1-5],  Hall  [2]^ 

Napar  stek  [1],  and  Petrovanu  [12]  for  even  more  complete  results 
as  well  as  additional  references. 

As  another  illustration,  we  briefly  describe  the  ideas  in 
the  paper  of  Landesman  and  Lazer  [1] .  Let  D  be  a  bounded  domain  in 
Rn  and  let 


n 


A  - 


i,  0=1 


a 

ck. 


10 


& 

dx". 

0 


be  a  second  order,  self-adjoint  uniformly  elliptic  operator  on  D* 
that  is,  each  a1^  =  a^1  is  assumed  to  be  real,  bounded  and  measurable 
on  D  and  there  is  a  constant  c  >  0  such  that 

n  n 

E  a10 1  |  >  c  E  £7 

i,  j=l  J  ”  0=1  J 


o 

for  all  x  e  D  and  for  all  real  £ . .  Let  H,  be  the  closure  through 

0 

the  (real)  inner  product 
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<u  v>  =  /  (u-v+u  -v  )dx 
1  D  X  X 

2 

of  C1  functions  with  compact  support  in  D;  let  Hq  =  L  (D)  where, 
for  f,g  e  Hq, 

<f,g>0  =  /  fgdx. 

D 

For  h  e  Hq,  a  >  0  constant,  and  g!  R  ->  R  bounded  and  continuous, 
the  problem  is  to  determine  a  weak  solution  of 


Au  +  au  +  g(u)  =  h(x),  x  e  D’ 
u(x)  =  0,  x  €  ciD. 


(7.12) 


That  is,  if 


n 


B(u,v)  =  Z  /  aljux  vx  dx, 
i,  j=l  D  i  j 


we  seek  a  u  e  such  that 


def 


(cp,Au)  =  B(cp,u)  =  <Jp,h-au-g(u)> 


O' 


(7-13) 


o 

for  all  cp  e  H-j- .  For  any  u  e  Hq,  let  Hu  =  h  -  g(u) .  Equation 
(7.13)  will  be  written  symbolically  as 


(A+al)u  =  Nu. 


(7*14) 
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With  the  representation  (7.13)  for  the  operator  A,  we  can 
o 

consider  A'  2>{A)  -  ^CHq  -»Hq.  A  special  case  of  the  result  of 
Landesman  and  Lazer  [1]  is  the  following! 

Theorem  7.3.  Suppose  9t(A+oil)  is  one-dimensional  and  w  is  a 
basis  for  ^(A+al),  w(x)  >  0  for  x  e  D,  Iw1q  =  !•  Assume  also 
that  g(s)  -»  g(«),  g(-s)  -»  g(-»)  as  s  ->  «,  g(«>),  g(-«)  finite, 
and 

g(-oo)  <  g(s)  <  g(»),  s  €  (-“»,«).  (7.15) 

A  necessary  and  sufficient  condition  for  the  existence  of  a  weak 
solution  of  (7-12)  is 

g(-«)  <  <h,w>Q  <  g(»).  (7.16) 

Sketch  of  Proof!  The  Riesz  representation  theorem  and  Rellich's 
compactness  theorem  imply  that  &(A)  =  Hq  and  that  A  has  a 
completely  continuous  right  inverse  S!  Hq  -»Hq.  Equation  (7-l4)  is 
equivalent  to  the  equation 

u  =  -aSu  -  S[g(u)  -  h]. 

Since  S  is  symmetric,  one  easily  obtains  <w, g(u)  -  h>Q  =  0,  and, 
from  the  fact  that  w(x)  >  0  in  D,  this  implies 
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g(-oo)  <  <h,w>Q  <  g(») 


is  necessarily  satisfied. 

The  proof  of  sufficiency  is  much  more  complicated  and  uses 
the  ideas  mentioned  in  Section  2.  Let  ^(A+al)  =  where  Q  is 

the  orthogonal  projection  through  <  • , • >^  onto  the  suhspace  spanned 
by  w.  The  operator  A  +  al  has  a  completely  continuous  right  in¬ 
verse  Ml  Xq  j  ^  ->Xq  j  Equation  (7.14)  is,  therefore,  equivalent 
to  the  equations 


u  =  Qu  +  M(l-Q)Nu; 
QNu  =  0. 


(7-17) 


If  we  let  Qu  =  aw,  where  a  is  a  scalar,  then  these  equations  can 
be  written  as 


P  ~P 

u  =  aw  +  M(l-Q)Nu  =  G-.  (u,a) 

(7-18) 

dpf 

a  =  a  -  <Mu,w>0  =  G2(u,a). 

If  we  let  G  =  (G1,G2),  then  finding  a  fixed  point  of  G  in  HQ  X  R 
is  equivalent  to  solving  the  problem.  The  difficult  part  of  the 
proof  remains.  By  using  the  hypotheses  on  g  and  some  careful 
estimates,  the  authors  show  by  means  of  Schauder's  theorem  that  a 
fixed  point  of  G  exists. 

Using  techniques  more  closely  related  to  those  of  Lazer  and 


-62- 


Sanchez  [1]  discussed  in  Section  4,  Landesman  and  Lazer  have 
also  proved  the  following  interesting  result!  Suppose  D  is  a 
domain  whose  boundary  is  regular  enough  for  the  existence  of  the 
Green' s  function  for  the  boundary  value  problem 

Au  =  h(x),  x  e  D,  u(x)  =  0,  x  e  cCD' 

where  A  is  the  n-dimensional  Laplacian  and  h  is  Holder  continuous 
on  D.  Suppose  that  for  <  a  <  there  exists  no  nontrivial 
solution  of 


Aq  +  au  =  0,  x  €  Dj  u(x)  -  0,  x  e  c5D. 

It  is  then  concluded  that  for  F  a  real-valued  C-^  function  on  R 
satisfying  <  F1 (u)  <  a g  for  all  u,  and  for  h  Holder  continuous 
on  D  U  3D,  the  boundary  value  problem 

Au  +  F(u)  =  h(x),  x  e  D*  u(x)  =  0,  x  e  £d] 

has  a  unique  solution. 

For  other  results  on  elliptic  boundary  value  problems  using 
the  ideas  of  Section  2,  see  Cesari  [10]. 
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8.  Admissibility.  In  this  section,  we  indicate  in  a  few  remarks  how 
the  important  concept  of  admissibility  as  defined  by  Antosiewicz  [2] 
subsumes  the  ideas  of  Section  2.  Suppose  X, Y  are  Banach  spaces,  F 
is  a  vector  space  with  at  least  the  structure  of  a  complete  topological 
vector  space,  G  is  a  subspace  of  F  with  at  least  the  structure  of 
a  Banach  space  with  the  topology  of  G  being  finer  than  that  induced 
by  F.  Let  u!  X  -» F  be  a  homomorphism,  v!  E  -» F  is  a  continuous, 
linear  transformation,  w!  G  ->  Y  is  a  continuous  transformation.  The 
first  problem  is  to  find  an  x  e  X,  g  e  G,  such  that 

g  =  u(x)  +  v  °w(g) .  (8.1) 

We  say  (Y,G)  is  admissible  (relative  to  the  pair  of  linear 
transformation  (u,v)  if,  for  any  y  €  Y,  there  is  an  x  e  X  such 
that  u(x)  +  v(y)  e  G.  An  application  of  the  closed  graph  theorem 
gives 

Lemma  8.1.  If  (Y,G)  is  admissible,  there  exists  a  p  >  0  and  for 
each  y  e  Y  a  point  x  e  X  such  that  u(x)  +  v(y)  e  G  and 

p||u(x)  +  v(y)||  <  || y|| . 

Lemma  8.1  is  basic  for  the  discussion  of  the  existence  of 
solutions  of  the  nonlinear  problem  (8.1)  when  w  is  small  in  some 
sense.  More  specifically,  for  any  given  g  e  G,  one  can  find  an 


-64- 


x  £  X  such  that  H(g)  =  u(x)  +  v°w(g)  e  G  and  HI  G  -»  G  will  he 

continuous .  If  one  can  find  a  subset  G  of  G  such  that  H*.  G  -» G 

o  o  o 

and  H  has  the  fixed  point  property^  then  there  will  be  a  solution  of 
(8.1)  in  G  .  A  result  of  this  type  is  given  by  Antosiewicz  [2], 
Applications  of  this  idea  to  the  theory  of  asymptotic  behavior  of 
solutions  of  ordinary  differential  equations  may  be  found  in  Hartman 
and  Onuchic  [1],  Implications  in  ordinary  differential  equations  of 
an  hypothesis  of  admissibility  on  the  behavior  of  the  solutions  of  a 
homogeneous  equation  may  be  found  in  Mas sera  [1],  Mas sera  and 
Schaffer  [1], 

It  is  easily  seen  that  the  problem  of  Section  2  is  a  special 
case  of  the  general  problem  considered  here.  Suppose  X,Z  are  Banach 
spaces  A,N  as  in  Section  2}  9t(A)  =  Xy,  j^(A)  =  Z_.  Then  Ax  =  Nx 

is  equivalent  to  the  system 


a)  x  =  Ux  +  MENx 

b)  (i-E)Nx  =  0. 


(8.2) 


If  Y  =^(A),  X  =  G  =  F,  u  =  U,  v  =  M,  w  =  EN,  then  the  set  (Y,X) 
is  admissible  by  construction  and  the  first  equation  (8.2a)  is  a  special 
case  of  (8.1) . 

The  concept  of  admissibility  is  a  sweeping  generalization  of 
the  Fredholm  alternative  and  will  continue  to  have  significant  applica¬ 
tions  to  differential  equations . 
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